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^ ■ 1. Introduction 

\ In this paper we solve the nonlinear differential equation 

^ : Yit? = (1 1) 

o: 

CSJ . where f{x) = cqx^ + 4aix^ + 6a2X^ + Aa^x + 04 is a quartic polynomial with no repeated 

\ factors and n > is a fixed whole number. The solution is expressed parametrically in 

^ ■ terms of the Weierstrass p-function p{w), and his sigma and zeta functions a{w), C{w); see 

formulas (2.4), (2.5), (2.6), (2.7), and definitions (A.l), (A. 5) of the Appendix. One can 
also solve equation (1.1) in case f{x) does have a repeated factor - a situation which is 
■ - - ■ easier to deal with and which, in particular, generally does not involve elliptic functions; 

see Example 6. 

Of some special interest is the equation 

Y{tf = BYitf + EY{t) -K + ^^+ (1.2) 

with n = 1 in (1.1). Note that the solution of the differential equation 

Yi{tf = Yi{tf [Bi + EiYiHr - KiYiitf"" + AiYiitf"" + DiYiit)''"'] , (1.3) 

where m 7^ is a fixed whole number, can be obtained from that of p.2|) . Namely, if we set 
Y{t) = yi(t)"™ then equation is transformed to ((121) for B = m?Bi,E = m?Ei,K = 
m'^Ki,A = m'^Ai, and D = m?Di. 



1 



21,2012 1:11 WSPC/INSTRUCTION FILE ParametricWS 



2 Jennie D'Ambrotse and Floyd L. Williams 

If both E and D are zero, for example, then the general solution Ye.d=o of ()1.2p . in 
terms of the Weierstrass elliptic function p{w), was obtained in 1933 by G. Lemaitre [18j . 
in his study of spherically symmetric distributions of matter. Compare also the paper [2D] 
of G. Omer with references therein to special case solutions by R. Tolman, B. Datta, and 
H. Bondi. The solution Ye,d=o provides for an exact, inhomogeneous cosmological solution 
ds^ of the Einstein field equations with cosmological constant A = ZB. Namely, for the 
family of Szekeres-Szafron solutions 



the functions B{x,y, z,t), A{x,y, z,t) are explicated by the solution Ye^d=o- The reference 
|15| . for example, contains a detailed discussion of this matter, and the text [16] can be 
consulted for a comprehensive analysis of inhomogeneous cosmology. Also see [T7], where 
one can consider D = -Q^ ^ q (with yet ^ = 0), Q being a constant electric charge. 

The Friedmann-Lemaitre-Robertson- Walker (FLRW) metric can be obtained from ds'^ 
in ()1.4p (as a "limit"), and a known formula (see formula (42) of [15], for example) for the 
scale factor (the "radius" of the FLRW universe) also follows from the general formulas 
presented here. Further remarks on this, as well as the computation of scale factors in 
anisotropic models (in the Bianchi V and IX models, for example) are taken up in section 



Another case of interest, among others to be mentioned later, is that when both A 
and B are zero in ()1.2p . Here the solution Ya,b=o is of relevance regarding the dynamic 
correspondence between Bose-Einstein condensates (BECs) and FLRW/Bianchi / cosmol- 
ogy |1HI19]. In particular we deduce an alternate formula (in section 3) for the second 
moment l2{t) = ^,B=o(i) of the wavefunction of the Gross-Pitaevskii equation, for BECs 
governed by a time-dependent, harmonic trapping potential - especially when a cosmological 
constant is present - i.e. when E ^ 0. The second moment in fact determines the harmonic 
trapping frequency. 

Our formulas therefore provide for a general, unifying context where some known formu- 
las in the literature are immediately derived or extended, and some new ones are developed 
- as seen in the concrete examples of section 3. 

2. The Solution Of Equation ([TIT]) 

For Z '^=' \J a{)Y'^ + ^a\Y'^ + 602^^ + ^a^^ + 04 equation (1.1) is expressed as 



for an integration constant 5. By a choice of any root xq of /(xq) = and a (finite) Taylor 
expansion /(x) = 4a3(x — xq) + 602 (x — xq)^ + 4ai(x — xq)^ + ao(x — xq)"^ of /(x) about xq, 
where as = /'(xo)/4, a2 = /"(2;o)/12, a\ = /"'(xo)/24, qq = /""(xo)/24 one can reduce 
the elliptic integral / in (2.1) to a Weierstrass canonical form, where the quartic in Z is 



ds^ = dt" - e2^(^'^'^'*)(dx2 + dy^) - e2^("'^'^'*)dz2 



(1.4) 



3. 




(2.1) 
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reduced to a cubic. Namely, by the substitution Y = xq + ol^/{x — 02/2) one gets 



/ = / dx (2.2) 



where 



g2 = — 4aia3 = 0904 — 4aia3 + Sa^, (2.3) 

(73 = 2010203 — O2 — OoQ^S = 000204 + 2010203 — O2 — O0O3 — 0]^04 

are the Weierstrass invariants of f{x). Here 03 7^ since f'{x) 7^ 0, as f{x) has no 
repeated factors, by hypothesis. Moreover if p{w) = p{w; g2, gs) is the Weierstrass p- 
function attached to §2,93 (as in the Appendix) then, by equation (A. 2) there, the substi- 
tution X = p(w + c), for any fixed constant c, leads to the re-statement 

xo + TT-r {^"""i, dw + 5 (2.4) 

A[p{w + c)- f"{x^)/2A]) ^ > 

of equation (2.1), where we also have (by the above substitution Y ^ x) 

'^° + 4[p(u; + c)-/"(xo)/24]- ^"''^ 

Equations (2.4), (2.5), which constitute the main result, provide for a parametric solution 
of equation (1.1). For the applications we have in mind, in examples 1-6 below, we need 
only the cases n = 0, 1, 2. For n = and the choice 5 = 0, t = w in (2.4) and formula (2.5) 
also follows from a general formula of Biermann- Weierstrass [6j; more on this point later. 
Let a{w)X{w) denote the Weierstrass sigma and zeta functions respectively; see definition 
(A. 5). For a choice wq with p{wq) = /"(xo)/24 we have for n = 1,2 respectively in (2.4) 



t = XqW + 



o-(u; + c + luo) 



+ <5, (2.6) 



x\w + 



x^r{x^) , /'(xo)V"(^«o) 



+ 



2p'{wq) 16p'(wo)' 



log 



a{w + c + tDo) 



o{w + C — Wq) 



16p'(wo)^ 
+{w + c 



[C(w + C + ZDo) + C(^^ + C - Wq) 



(2.7) 



a:o/ (xq) ^^ /(xq)^ 
— 77 — T-C(''«o) 77 — 



2p(u;o) ^ 2p"iwo)C{wo) 



p'(u;o)^ 



+ 5, 



by formulas 1037.06, 1037.11, respectively, in [7j - assuming also that p{wo) 7^ the roots 
ei, 62, 63 of 4x3 — 523; — 93 = 0. In fact, in principle, one can compute the integral in (2.4) 
for arbitrary n > 0, using formulas in |7j, which we shall not pursue as the result is not 
needed here. 
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The aforementioned Biermann-Weierstrass formula for the solution Y{t) of equation 
(1.1) in case n = is 

/(y(o))VV'(t) + m (pit) - mi) + fiymrivm - 
Y{t) = Y{Q) + ^ , 1^, {, (2. 



P{t) - ^ 



/(y(o))/""(y(o)) 

48 



Also see references [2TII22]. 



3. Some examples 

The following examples are meant to provide application and further clarity of the preceding 
formulas. 

Example 3.1. We begin with equation (20) 

■ . x9 .X 2M(z) A^(tf , , 

$(t)2 = -K{z) + + (3.1) 

of |15j . which is associated with the metric ()1.4p in the introduction. Since the functions 
K{z)^M{z) are independent of t, equation ()3.ip is mathematically the same as equation 
(8.9) of [18j, and since its solution is presented and discussed in [15] we greatly limit our 
remarks here. (jS.ip is equation ()1.2p with 35 = A (a cosmological constant), E = 0, K = 
K{z), A = 2M{z),D = 0. Thus /(x) = (A/3) - K{z)x^ + 2M{z)x has xq = as a first 
order root for M{z) + 0, since /'(O) = 2M{z). By formulas (2.5), (2.6) the solution of (i3J]l 
is given parametrically by 



W^ + C;ff2,<?3)+i^(-2)/12 

(3.2) 



M{z) 
2p'{wo) 



a(w + c — Wo) , X X 
o-(u; + c + -Wo) 



+ 5 



for p(u)o) = — -ftr(2;)/12 / 61,62,63, which are formulas (21), (23) of [15], where g2 = 
K{zf/12,g3 = K{zf / 216 -AM (z)^/ 12 (by (2.3)), and where now the integration constant 
6 = 5{z) depends on the variable z. 

Example 3.2. Consider the Friedmann equation 

a'ivf = ->^a{vf + ^ + ^a(r,) + ^a(ry)^ (3.3) 

3 

for the scale factor a(r/) in a FLRW universe whose energy and matter are modeled by 
a perfect fluid. Here rj is conformal time, k is the curvature parameter, K = SttG/S for 
the Newton constant G, Aj. and A^ are radiation and matter constants, c is the speed of 
light, and A is a cosmological constant. Given the "big bang" initial condition a(0) = one 
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derives immediately from (2.8) the result 

^ 2 (p(r/) + k/12)2 - KAAr/Qc^ ' 

with known application in the study of cosmic microwave background fluctuation, for exam- 
ple; see [2l[3]. The invariants §2,93 in (2-3) are given by 



(3.5) 



3c2 + 12 



~ 18c2 "^ 216 48c4 
Elliptic function solutions of Friedmann equations are also discussed in p!1[9lll3j . for example. 

Example 3.3. As another example we consider the Bianchi V cosmological model with 
metric 

ds^ = -dt^ + X{tfdx^ + e^'"'Y{tfdy'^ + e^^"" Z {tf dz'^ (3.6) 

for 6 7^ 0. As before, we take the energy momentum tensor to be that of a perfect fluid 
and denote the radiation and matter constants by Ay^ Am, respectively. For K = SttG/S, a 
zero cosmological constant A = 0, and the speed of light taken to be c = 1, the Einstein 
equation is a special case of (1.1) 

Ritf = [b^Rit)^ + KAmR{tf + KArR{tf + KD] (3.7) 

in terms of R{t) {X{t)Y{t)Z{t))^/^ and the quantity -D ' ^ (|| + ^ + || - 



— — Yzj '^tiich can be shown to be a constant; see [10]. Then by (2.5) and (2.7) we 
obtain 



4 \p{w + c) - /"(xo)/24] 

(3.8) 

2 , ( xof'ixo) , f'{xo)'^p"{wo)\ a{w + c + wo) 

: ^ : loff : 



* ''o^+y 2p'(t«o) ^ IQp'iwof )^''^a{w + c-wo) 

^^"'^ ^2 ['^("^ + C + Wo)+ C(.W + C - Wo)] 



..r. f'(..r.\ f'(..)^ r 2p(u;o) ^ 2p"iwo)Ciwo) ] \ ^ ^ 

H 77 — +d 



, I , X / Xpf (xq) /'(xo) 



P'{wq) 16 



P'{wqY p'{wof 
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for p{wq) /"(a:o)/24 and associated polynomial f{x) = b'^x^ + KA^x^ + KArX^ + KD. 
The invariants are 

(3.9) 

,3 = :^f,2^,0 ^'rX 3D/1?„A- 



6 V 36 8 

where as in the Appendix we assume that — 27(7| ^ 0. To obtain the metric one writes 
X = R,Y = Re^^^,Z = Re-^^^ where r(t) = / dt/R{tf. That is, r is given 
parametrically in terms of w by 



w f'{xo) 



aiw + c — wi) , . ^, s 



+ 5' (3.10) 



for integration constant 5' , p{wi) = f"{xo)/2 — f'{xo)/Axo and assmning xq ^ 0. 

Example 3.4. For a cosmological model with Bianchi IX metric 

ds^ = -dt'^ + a{tfdx'^ + b{tfdy'^+[b{tfsm'^y + a{t)'^cos^y] dz^ -2a{tf cosy dxdz, (3.11) 

and massless scalar field <I> and flat (constant) potential V{r) = 2A, the modified Einstein 
equations based on Lyra geometry are studied in [3], for example. Also compare the paper 
[3]. The field equations yield the relation ^ = <I>o/a6^ for an integration constant $0) and 
the assumption a = leads to the Einstein equation 

b , .b"^ 1 62n-4 

7+n + lW = - — 3.12 

b b"^ [n — 1)6^ n — I 

for n 7^ 1, which moreover is shown to have the first integral 

for n ^ 0, ±1, where Di is an integration constant. The authors found solutions of ()3.13p 
only for Di = - for n = 2,1/2,3/2,3/4. Therefore we consider the case Di ^ 0, and we 
take n = 2 for example: 6^ = 1/3 - 6^/4 + Di/b^, or a = 5^ 

a{tf = ^a{t)-a{tf + ^, (3.14) 

which is another example of equation ()1.2p for D = 4Di,A = 0,K = 0,E = = — 1. 
Here 

/(x) = -x^ + ^x3 + 4Z)i, (3.15) 
g2 = -4L>i,53 = — = — . 

In particular ~ 27(73 ~ —16Di{4:Di + 1) ^ for Di ^ 0, — which we assume. Certainly 
if Di = — then /(x) = — x'^ + |x^ — | = — (x — 1)^ (x^ + |x + ^) has xq = 1 as a repeated 
root, for example. 
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On the other hand, consider xq = —1 which is a non-repeated root for Di = j2 '■ 
/'(xo) = 8/0. By (2.5), (2.6) we have the parametric solution 

9 

a = -1 + 



t = -w + 



p{w + c) + 5/6' 

(3.16) 

2 " 



p'iwo) 



a(w+C — WQ) ^ 

log -j ■ + 2{w + c)Ciwo) 

a[w + c + wo) 



+ 6 



of equation ()3.14p for p{wq) = —5/6. One can find solutions to equation ()3.14p for many 
other non-zero values of Di similarly. This amounts to specifying Di 7^ —j2 and solving for 
a corresponding root xq of f{x) in ()3.15p . 

Example 3.5. The differential equation 

Y{tf = EY{t)-K+^, (3.17) 

which is another example of equation (jl.ip (for n = 1), arises in the study of time-dependent, 
harmonically trapped Bose-Einstein condensates, as indicated in the introduction. The con- 
stant E here is a positive multiple of a d— dimensional cosmological constant, for d > 3 
arbitrary. Also D is positive and E ^ 0. See equation (35) in frii . where some elliptic 
function solutions are discussed. Equations (2.5), (2.6) provide for a parametric solution 
of ()3.17p . However, as f{x) = Ex'^ — Kx^ + D is cubic in this case, one has an alternate, 
simpler parametrization which in particular does not involve a logarithm, as in (2.6). 

Namely consider the simple substitution Y = ax + b, a ^ 0, a, suggestion for which we 
thank one of the referees. Equation (2.1) then reads 

(a^x + ab)dx 

^Ea^x^ + {ZEaH - Ka^)x'^ + {ZEab'^ - 2Kab)x + Eb^ - KP + D 

(3.18) 

1 /■ (u + ab)du 



i2 



V^4m3 - g2U- 93 



for u = a^x,a = {E/AY/'^,b = K/3E,g'^ = {-3Eb^ + 2Kb)/a, and g'^ = -Eb^ + Kb^ - D. 
Note that 52 '53 the invariants 52 = -f^^/12,53 = ^ — ^f^ of /(x) given in definition 

(2.3). Similar to the derivation of equation (2.4), we now let u = p{w) = p{w; g'2, g'^) to get 
the parametric solution 

t = ^C{w) + -w + 5 



(3.19) 

of equation (|3.17p . where again C,{w) is the Weierstrass zeta function of (A. 5). 



Y = -p{w) + b 
a 



Example 3.6. In the main result, equations (2.4), (2.5), f{x) is assumed to have no 
repeated factors. However, if repeated factors occur then equation (1.1) can be solved, 
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in fact, in terms of elementary, non-elliptic functions, which we illustrate in the following 
example. 

The equation 

U'ixf + Wixf - 2U{xf - U(x)/V2 = — , (3.20) 

16 

an example of equation (1.1) with n = 0, is satisfied by the potential U{x) of a 
particular Zakharov-Shabat system. U{x), moreover, satisfies a type of static, modified 
Novikov-Veselov (mNV) equation [14|, a point which we return to later. Here f{x) = 

— + 2x^ + x/V2 + jQ has (x — xq) as a repeated factor for xq '^=' — -v/2/4, f{x) = 
(x — xo)^ [Bi + B2(x — Xq) + B^^x — xq)^] being its finite Taylor expansion about xq for 
2Bi = /"(xo),6B2 = /'"(xo),2453 = /""(xq). By the substitution x = xo + ^, 

I'M- [^ = ±f ^ (3.21) 

is an elementary integral (compare the integral in (2.1)) whose evaluation depends on the 

signs of Bi and the discriminant A '^=' B2— ^BiB^. Actually, as xq = — \/2/4, Bi = —1 < 
and A = 16 > 0. 
In the end, for 

1 . ( B2 + 2B^u \ 
I = arcsm — (3.22) 

and an integration constant 6 we obtain the solution 

Uix) = U^{x;6)= "^^Mx - 5) ^ ^ 

^ ' ' 2V2[V2±sin{x-5)] ^ ' 

of equation (3.20). In particular, one has the solution C/+(x;0) obtained in [H] (by a quite 
different method) where among other results the authors there demonstrate invariance of 
the Willmore functional W (the Polyakov extrinsic string action) under NV deformations. 
We note that the "deformation" 

, . def. def. sinjx + 2t) 

Uix,t) = U,ix + 2t;0) = ^7f|7|r— (^T^ (3.24) 

of f/+(x; 0) is a solution of the mNV equation 

Ut = U^^^ + 24U^U^. (3.25) 
The functional W of course is a basic quantity in the study of two-dimensional gravity. 



Appendix A. 

Given the central importance of the Weierstrass p— function p{w) for the present work we 
recall briefly, for the reader's convenience, its construction/definition. As we have indicated 
in section 2 a more detailed account is available in 

Let uji,uj2 be non-zero complex numbers. Since the imaginary parts of a non-zero com- 
plex number z and its reciprocal are related by Im = — (Im z) \z\~'^, one has that 
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Im UJ2/(^i 7^ if and only if Im uji/uj2 0. In particular we assume that Im uj^I'jJx > 0, 
which is equivalent to the assumption Im uji/uj2 < 0. The corresponding lattice ^ = 
^{(jji,uj2) generated by uji and UJ2 is defined to be the set of points uj = muji + nu;2 where 
m and n vary over the set of whole numbers. The lattice ^ gives rise to the p— function 



1 1 



(A.l) 



which is also denoted by p(w\ ), or by p{w] wi, a;2). p{w) is a meromorphic function, which 
is doubly periodic with periods wi, W2- Thus, by definition, p[vS) is an elliptic function. p{w) 
has double poles at w; = a; £ =Sf , and it satisfies the differential equation 

p'{wf =4p{wf - g2{u;i,uj2)p{w) -53(^1,^^2) (A. 2) 

for invariants 



<72K,^2) =^-60 Yl 53(u;i,u;2) =^-140 -1 (A.3) 

where, moreover. 



uje^-{0} iueJf-{o} 



52(^1,^2)^ - 2753(^1,^2) V 0. (A.4) 

Conversely, it is an amazing fact that if two complex numbers g2 and are given that 
satisfy the condition i^l — 27^(1 7^ 0, then there exists a pair of non-zero complex numbers 
iOi,uj2 with Im UJ2/UJ1 > such that 5'2('^i5'^2) = 92 and g3{uJi,UJ2) = 93, for 52(^1, (^2) and 
(73(a;i,cj2) defined in (A.3) with respect to the lattice = ^{uji,uj2) generated by uji and 
L02- Thus from g2 and g^ one can also construct the corresponding p— function p{'w; UJi , UJ2) 
(according to definition (A.l)), which in this case we also denote by p{w; g2, g-i) - as we 
have so done in the previous sections. 

Associated with p{w) are the Weierstrass sigma and zeta functions a{w) and C(^)5 
respectively: 



C'{w) -piw), hm (Ciw) - -) 0, 



a'{w) def. cr{w) def. 

— C[w), hm = 1. 



(A.5) 



a{w) ' w^Q w 
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